Spin-flip scattering at quantum Hall transition 
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We formulate a generalized Chalker-Coddington network model that describes the effect of nuclear 
spins on the two-dimensional electron gas in the quantum Hall regime. We find exact analytical 
expression for the transmission coefficient of a charged particle through a saddle point potential 
in presence of perpendicular magnetic field that takes into account spin-flip processes. Spin-flip 
scattering creates a metallic state in a finite range around the critical energy of quantum Hall 
transition. As a result we find that the usual insulating phases with Hall conductance o~ xy = 0, 1, 2 
are separated by novel metallic phases. 



Transition between the quantized values of Hall con- 
ductance in integer quantum Hall effect represents a pro- 
found example of electronic delocalization. Properties 
of this transition remain a subject of intensive theoret- 
ical and experimental investigations since the discovery 
of the integer quantum Hall effect (QHE) [lj]. Theory 
of the delocalization transition in the QHE predicts the 
existence of a single delocalized critical electronic state 
at the center of Landau level (LL) Q . Other states of a 
given LL are localized by disorder due to the Anderson 
localization phenomenon. Effective theoretical treatment 
of QHE delocalization transition is provided in frame 
of the Chalker-Coddington (CC) network model [3j. In 
this model electrons move along unidirectional links that 
form closed loops in analogy with semiclassical motion on 
contours of constant potential. Scattering between links 
is allowed at nodes, in analogy with tunneling through 
saddle-point potentials in the semiclassical model. In 
Ref. [4|] B. Shapiro pointed out that a convenient mod- 
eling of disordered systems can be given by networks of 
unitary random scattering matrices that correspond to 
the symmetries of the Hamiltonian. In this way, subse- 
quent generalizations of CC model allowed to include spin 
degree of freedom [j| and to describe systems belonging 
to various symmetry classes 0, 0] • 

In this paper we study the effect of spin-flip scattering 
by magnetic nuclei on the QHE transition. The spin-flip 
scattering by nuclei mixes two components of electron 
spin, and results in the mutual influence of electron and 
nuclear spins. The interplay of electron and nuclear spin 
degrees of freedom has been investigated earlier, focusing 
on measurement of nuclear-spin-lattice relaxation medi- 
ated by hyperfine interaction with two-dimensional elec- 
tron system in QH regime In contradistinction, we 
concentrate on the effect of nuclear spins on the electron 
motion and resulting changes in the QH phase diagram. 

We adopt the model of point-like exchange interaction 
between nuclear and electron spins H mt = JI • s, where I 
and s denote the spins of the nucleus and of the electron 
respectively. Throughout the paper we assume spin-1/2 
nuclei. In the absence of spin-flip scattering there are 



two Zeeman-split critical energies for each Landau level, 
where the QH delocalization transition takes place. We 
find that the spin-flip scattering results in the appearance 
of a finite region of delocalized states around the critical 
QHE states. The energetical width of the delocalized 
region depends on the strength of the exchange interac- 
tion relative to the Zeeman splitting in a nonmonotonous 
way. It increases with the exchange strength for weak 
exchange, reaches its maximum, and then shrinks back 
to two critical states at very large exchange interaction. 
Our results are summarized in phase diagram shown in 
Fig. [1] To obtain this phase diagram, we constructed an 
effective network model based on the famous CC model 
for the QH transition. 

Exact description of quantum Hall transition that 
takes into account scattering of electrons by impurity 
spins requires the knowledge of a many particle wave 
function ^(t, r, R4, ...Rjy), where r denotes the coordi- 
nate of the mobile electron and Rj , z = 1 , . . . , N are the 
coordinates of nuclear spins. It follows that a formal CC 
model has to use scattering matrices for the many parti- 
cle wave function introduced above. However, the short 
range of interaction between the electron and a nuclear 
spin allows for a series of simplifications that lead to an 
effective description in frame of CC model formulated for 
the reduced two-particle wave function. First, let us no- 
tice that far from the saddle points of disorder potential, 
the energies of spin-up and spin-down electrons are dif- 
ferent due to Zeeman splitting, and the elastic spin-flip 
scattering is suppressed. The energy conservation can 
however be respected near the saddle-points [l(| • There- 
fore, in terms of CC model, the spin-flip scattering is 
effective only at the nodes of the network. Second, the 
act of scattering in each node involves only the coordi- 
nates of the electron and a single nucleus that is located 
in the node. Considered alone, this scattering event can 
be described in terms of the scattering matrix for the two 
particle wave function *&i(t, r, R^), where r is the coor- 
dinate of the electron and Ri is the coordinate of the 
nucleus in the node i. 



v&° ut (^r,R 4 ) = ^ n (t,r,R 4 ) 
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FIG. 1: (Color online) a) Phase diagram at zero Zeeman split- 
ting, b) Phase diagram with Zeeman splitting. 



The explicit expression for the scattering matrix elements 
is obtained along the lines of Ref. [TO] generalized on 
the problem with spin-flip scattering. We assume, that 
the scattering at saddle points of the potential (that 
correspond to the nodes in CC model) is accompanied 
by the interaction with nuclear spin, and hence allows 
for spin flips. Thereby, the states | let/) and | idi) 
(where indices e,I correspond to the spin state of elec- 
tron and nucleus respectively) cannot flip the spin be- 
cause of conservation of the total angular momentum by 
scattering. The scattering matrix for these two states 



can be obtained directly from the expression of Ref. 1 1] 
by shifting the hight of the saddle point by Zeeman en- 
ergy of the electron. Here we neglect Zeeman energy of 
a nucleus, which does not change the results of this pa- 
per qualitatively. We note in passing, that the model 
that takes into account the Zeeman energy of the local- 
ized spin scatterers can also be applied to the scatter- 
ing by localized magnetic impurities as it takes place, 
for example, in semimagnetic semiconductors [lij . The 



states | Tel/) and | | e T/) can be mixed in the spin- 
flip process. To obtain the transmission coefficients for 
those states, we find the eigenstates and eigenenergies 
of the Hamiltonian H = H + H int + Hz, where H 
is the Hamiltonian describing the motion of the elec- 
tron in the scalar saddle point potential [ll|, and Hz 
describes the Zeeman energy. We choose the basis of sin- 
glet and triplet states that are the eigenstates of the in- 
teraction part of the Hamiltonian. The two-dimensional 
Hilbcrt space for the spin-flip scattering problem is then 
formed by the states * ,o = ^ (I Te, In) + | ie, Tat)) and 
*io = 75 (I U,In) - I IcTjv)). Similarly to Ref. [Il| 
we seek the solution of the Schrodinger equation in the 
form $„(X, s) = ^„(s)$„(X), where ip n (s) and $„(X) 
describe the cyclotron motion and the motion of the guid- 
ing center respectively, and n denotes the number of Lan- 
dau level (LL). Explicit expression for ipn( s ) is provided 
in Ref. . Without the loss of generality we restrict the 
further consideration to the lowest Landau level and sup- 
press the Landau level index. The equation for the guid- 
ing center wave function in the basis <f> = (*o,o, ^i,o) T 
is given by 
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We seek for the solution of Eq. ([2]) in the form 



$P0 = 0(x) 



(2) 
(3) 



where 4>{X) describes the spatial dependence of the wave 
function, and (Ci, C2) 2 " is its spinor part. Let 4>{X) be the 
solution of the Schrodinger equation 



H 4>{X) = E4>{X). 



(4) 



Diagonalization of the spinor part of Eq. ^/ for a fixed 
E leads to the two eigenvectors and corresponding eigen- 
ergies 
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£1,2 = E-^J±^^P 
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The energies Si correspond to the total energies of the 
two-particle states. By the energy conservation, both 
energies are to be equal to the energy of the incoming 
electron, E\ — £2 — E. This leads to the two differ- 
ent values for the energy E that determines the effective 
height of the potential barrier 
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Following [lll | , we introduce the dimensionless measure of 
energy e = {E+\ J) / Ei , where E\ is the energetic param- 
eter characterizing the form of the saddle point poten- 
tial. Furthermore, we denote the dimensionless strength 
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of the Zeeman coupling as A = B jE\ , and the relative 
exchange strength as S — J/E\. Therefore, there are two 
solutions of Eq. @ corresponding to the energy E of 
the incoming electron, \$ lj2 (X)) = 0(ei,2, X) |spin 1 2 ), 
where £1,2 = e T (1/2)VA 2 + <5 2 , and (f>(e lt2 ,X) is the 
solution of Eq. for the energy E lt2 respectively. The 
explicit form of (f> n (ti,2, X) is given in Ref. (Tlj. The spin 
parts of the wave functions are given by 
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where V = y/A 2 + S 2 . The energies ei,2 together with 
the energies e^i = e ± A/2 — (5/2 of the states |Tet/) 
and ||e|j) determine the transmission coefficients t(ti) 
defined below. From Eq. (JT]), we obtain the transmission 
coefficients for the scattering with and without spin flip: 

° Ut (ieU I Uiv) in = OUt (UlN I | e TiV> in 

= ±[t( ei )-t(e 2 )}, (8) 

V + A . , 2? -A . , 

+ ^p-*( £ 2)> ( 9 ) 

^.) + ^), do) 
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where from Ref. [11 1 



y/T+ cxp(-7rei) 
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The elements for reflection are given by the exchange 
t(a) — > r(ej) = \/l — t 2 (fi) everywhere. Using the ex- 
pressions for transmission and reflection coefficients, we 
construct the scattering matrix in the node relating the 
incoming and outgoing waves (see Fig. [2]). From that 
we derive the transfer matrix relating the waves on the 
left and on the right sides of the node. The system is, 
on average, invariant under 90° rotation if at the next 
neighbor node the transmission and reflection (of each 
channel) are interchanged, i.e. e$ —* — Q. 

In the standard formulation of the CC network model, 
the outgoing wave after the scattering at the node i and 
acquiring a random phase on the link plays the role of 
the incoming wave for the scattering by the next node 
(let's number it with i + 1). The latter allows to get the 
information about the propagation through the network 
by dividing the network into slices and multiplying the 
transfer matrices slice by slice. Such a direct approach 
is inapplicable if one replaces the exact many-particle 
scattering matrix by the two-particle scattering matrix 
though. Namely, the outgoing two-particle wave func- 
tion from the node i contains no information about the 
phase and the spin direction of the nucleus in the node 
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FIG. 2: (Color online) Incoming and outgoing states at a 
single node. 



z + l. This is in contradiction with the direct multipli- 
cation of the transfer matrices that would automatically 
make the outgoing state of the nucleus i equal to the in- 
coming state of the nucleus i + 1. Therefore, in order to 
use the two-particle scattering matrix, one has to mimic 
the loss of the information about the state of the nucleus 
i before the scattering by the nucleus i + 1. For exam- 
ple, the outgoing two-particle amplitude for the spin-up 
electron contains two contributions, | tet/) an d | Tell)- 
In order to erase the memory about the state of the pre- 
vious nuclear spin, we allow those two components to be 
exchanged randomly before being scattered by the next 
node. Similarly we exchange the components | let/) and 
I ieii). Effectively, we introduce a x Pauli matrix act- 
ing randomly with probability 1/2 on each link on the 
state of the nucleus. This mechanism can also be inter- 
preted as a spin-orbit scattering acting on nuclear spins 
only. Therefore, we arrive at the version of CC model 
that describes the propagation of a bi-spinor consisting 
of an electron and a nuclear spin through the network. 
In addition to a usual random phase contribution at each 
link, there is a spin-orbit scattering on the links that acts 
only on the nuclear component of the bi-spinor. In Refs. 
12L 1 1 31 ] it was shown that spin-orbit scattering can lead 
to the appearance of metallic phases in phase diagrams 
of network models. Our numerical simulations result in 
similar conclusions. As the Fermi energy crosses the re- 
gion of delocalized states, there is a smooth change of the 
Hall conductivity by one unit between the two quantized 
values. 

Our numerical calculations proceed as follows. We 
study a system of size M x L where M is the width 



4 



of the system and L is its length. For a given M, the 
eigenvalues of ln(T + T), where T is the full transfer ma- 
trix, behave as exp(— 2X n L) defining the Lyapunov ex- 
ponents A„; the smallest positive one, Xm/2j defines the 
localization length £m/M. The M dependence of £,m/M 
identifies the phases: (i) a decreasing ratio corresponds 
to localized state, i.e., an insulator, (ii) a constant ra- 
tio corresponds to a critical state, and (iii) an increasing 
ratio corresponds to a metallic phase. 

The results of our numerical calculations lead to phase 
diagrams presented in Fig. [1] The phase diagram in 
Fig. QJi corresponds to the spin-degenerate case with zero 
Zeeman splitting. In the absence of spin-flip scattering 
(5 = 0) there is a transition between the two insulating 
phases with Hall conductances a xy — and a xy — 2 at 
energy e = 0. At finite S, a metallic phase appears in the 
finite range of energies, as described above. 

The phase diagram for finite Zeeman splitting A = I 
is shown in Fig. [TJd. Without the spin-flip scattering 
there are two critical energies e c = ±A/2 separating the 
insulating phases with a xy = 0,1,2. A finite spin-flip 
scattering induces the appearance of two metallic regions 
around each critical energy, in analogy with Fig. As 
the spin-flip scattering rate grows further, the size of each 
metallic region diminishes and finally collapses to single 
critical states at energies e w ±A/2 + 6/2. 

The appearance and the subsequent collapse of dclo- 
calized phase with the increase of the spin-flip scattering 
5 can be qualitatively explained by the detailed consider- 
ation of the mutual influence of the spin-flip scattering at 
nodes and mixing of the states on links. Elements of the 
scattering matrix Eqs. © - (JTUJ) describe the scattering 
of the states with antiparallel spins of electron and nu- 
cleus. Propagating along the links, those states mediate 
the indirect mixing of the states | let/) an d | ieil) by the 
mechanism described above. A typical sequence of tran- 
sitions leading to a mixing between | TeTi) and | | e |/) 
looks like | let/) -> I Tel/) on the link, | Tel/) -* | leT/) 
at the node, | | e |/) - * | lei/) on the link. This mixing 
acts similarly to the spin-orbit interaction for the electron 
wave function and, if effective, results in the appearance 
of finite regions of dclocalized stated around the critical 
energies e-f and ej_. The effect of the mixing is reduced 
with the separation between the energies e-f and ej. , and 
ei , €2 ■ For small values of 5 there is a finite interval of en- 
ergies e, in which the transmission amplitudes t(ei) and 
£(£2) are of the same order. In this interval, the effec- 
tive spin-orbit interaction results in the metallic regions 
around the critical energies corresponding to the states 
I TeT/) an d I lei/)- At large S, the large energy separa- 



tion between e\ and ti leads to the localization of at least 
one of the states corresponding to those two energies for 
any e. As a result, the mixing on the links necessarily 
involves a localized state and hence becomes ineffective. 
The states | TeT/) an( i I lei/) become two decoupled 
critical states, the metallic regions collapse (see Fig. [I}. 

In conclusion, we derived the exact solution of spin-flip 
scattering problem in quantum Hall regime, and used it 
to construct the network model that describes integer 
QH transition in presence of scattering by nuclear spins. 
With increase of the exchange coupling between electron 
and nuclear spins, we predict the appearance of metallic 
regions in the QH phase diagram that subsequently col- 
lapse to the critical lines. Those critical lines correspond 
to a standard QH transition. Our investigation can also 
shed light on the QH transition in presence of scattering 
by magnetic impurities. 
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